      SUBROUTINE SETUP
      RETURN
      END
C   HS64
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=5.D0*X(1)+5.D4/X(1)+2.D1*X(2)+7.2D4/X(2)+1.D1*X(3)+1.44D5/X(3)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)= 5. D0-5.D4/X(1)**2
      GRADF(2)= 2. D1-7.2D4/X(2)**2
      GRADF(3)= 1. D1-1.44D5/X(3)**2
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 )CRES(I+NH)=CRES(I+NH)+1      
      IF(I.GT.3)      GOTO 100
      GXI=(X(I)-1.D-5)
      RETURN
  100 CONTINUE
      GXI=1.D0-4.D0/X(1)-32.D0/X(2)-120.D0/X(3)
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADGI(NX)
       IF ( I .GT. 3 ) GOTO 100
       DO 50 J=1,3
         GRADGI(J)=0.D0
   50  CONTINUE
       GRADGI(I)=1.D0
       RETURN
  100  CONTINUE
       CGRES(I+NH)=CGRES(I+NH)+1
       GRADGI(1)=4.D0/X(1)**2
       GRADGI(2)=32.D0/X(2)**2
       GRADGI(3)=120.D0/X(3)**2
      RETURN
      END
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS64'/
      DATA (X(I),I=1,3)/3*1.D0/
      DATA N/3/ , NH/0/ , NG/4/
      DATA DEL0/1.00D0/ ,TAU0/1.0D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,4)
     F      /-1,0,0,1,1,1,1,2,1,1,3,1,-1,0,0/
      END
